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QCD-like anomalous contributions to the dimensions of gauge theory operators is 
proposed. This approach leads to a relation between scalar glueball masses and 
beta functions. Using this relation, the glueball mass spectra that emerge from 
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functions are found. Remarkably, some of them present an infrared fixed point at 
finite coupling. 
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1. Introduction 

The idea that string theory and non-abelian gauge theories are related has been 
proposed long ago [1]. More recently, an exact equivalence between string theory in 
ten dimensions and gauge theories in four dimensions was discovered [2, 3, 4]. This 
relation holds for string theory in AdS 5 x S 5 spacetime and SU (N) Yang Mills theory 
with large N, extended Af = 4 supersymmetry and conformal invariance. This is 
an example of the AdS/CFT correspondence which also includes other gauge/string 
dualities in different geometries and dimensions. 

In order to apply the idea of gauge/gravity dualities to describe strong inter- 
actions, it is necessary to break the conformal invariance. Presently, an exact dual 
description of QCD is not known. However, in the recent years, some important 
QCD properties have been reproduced from phenomenological models based on the 
AdS/CFT correspondence. Essentially, these so called AdS/QCD models consist in 
modifying the AdS geometry with the purpose of breaking the conformal invariance. 
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An infrared scale in the gauge theory was associated with a localization in the 
AdS space in Ref. [5]. This way a physical process with an infrared scale in four 
dimensions is mapped into a region of the AdS space. Using this idea, the correct 
high energy scaling of hadronic amplitudes for fixed angle scattering was found. This 
experimentally observed scaling was reproduced by QCD long before in [6, 7, 8]. 

Putting forward the idea of Ref. [5] to relate AdS/CFT and QCD, in Ref. [9, 10] 
the scalar glueball spectroscopy was studied using an AdS slice. The size of the slice 
is related to Aqcd- Considering boundary conditions in the AdS slice, normalizable 
modes for the scalar bulk fields dual to scalar glueballs with a discrete spectrum were 
found. The approach of using an AdS slice to investigate hadronic properties was 
then called AdS/QCD Hard-Wall model and applied to other particles (see e.g [11]). 

It was then realized that the Regge trajectories that can be obtained from the 
Hard- Wall model are not linear. Then, another AdS/QCD model was proposed 
to give linear trajectories, especially for vector mesons. This was done with the 
introduction of a dilaton background field that plays the role of a smooth cut-off in 
the AdS spacetime [12]. This is called Soft- Wall model and was also applied to study 
scalar particle properties [13, 14, 15, 16]. 

Another interesting approach to the holographic description of QCD appeared 
in [17, 18] (see also [19] for a review). In these models, dilaton potentials are related 
to QCD beta functions. This makes it possible to use these holographic models 
to obtain many QCD properties. In particular, they obtained glueball masses in 
agreement with lattice results. 

In the AdS/CFT correspondence, the masses of supergravity fields are related 
to the scaling dimension of the local gauge-invariant dual operators. Since the gauge 
theory is conformal, the beta function vanishes. In the original AdS/QCD models, 
one usually assumes that the relation between the masses of the bulk fields and 
the scaling dimension of the boundary operators is the same as in the AdS/CFT 
correspondence. 

The purpose of this work is to study the scalar glueball mass spectrum taking 
into account contributions from the anomalous dimension of the glueball operator 
in a Soft- Wall background. This anomalous dimension implies a modification of the 
mass of the dual five dimensional supergravity field. This modication implies in turn 
a correction on the four dimensional mass spectrum. We study the case of the scalar 
glueballs, that are dual to scalar bulk fields in the holographic AdS spacetime. 

We consider some possible phenomenological beta functions and find the corre- 
sponding glueball mass spectra. We also investigate the inverse problem: starting 
from a linear fit of lattice quadratic glueball masses, we identify possible correspond- 
ing beta functions. 

The paper is organized as follows. In section 2, we review lattice results, discuss 
the Soft- Wall model and calculate the anomalous dimension of the scalar glueball 
operator. In section 3, we obtain glueball spectra from three different beta functions 
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considered in the literature. Section 4 is devoted to constructing beta functions able 
to reproduce linear Regge trajectories of lattice scalar glueball masses. 



2. Scalar Glueballs 
2.1 Lattice results 

The glueball mass spectrum has been the subject of numerous studies in lattice QCD 
[20, 21, 22, 23], which have provided estimates for the groundstate and a few excited 
states as shown in Table 1 (for general reviews, see for instance [24, 25]). 





Ref. [20] 


Ref. [21] 


Ref. [22] 


Ref. [23] 


jPC 


N c = 3 


N c = 3, anisotropic lattice 


N c = 3 


N c ->■ oo 


0++ 


1.475(30)(65) 


1.730(50)(80) 
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1.58(11) 


1.48(07) 


++* 


2.755(70)(120) 


2.670(180)(130) 




2.75(35) 


2.83(22) 


Q++** 


3.370(100)(150) 










Q-| — |-*** 


3.990(210)(180) 











Table 1: Lattice scalar glueball mass spectra in GeV. The errors are shown in parenthesis 
and described in the text. 



In the second column, we show the results of Ref. [20] computed for pure SU(3) C 
lattice gauge theory. The first number in parenthesis is the statistical error stem- 
ming from the continuum-limit extrapolation while the second error accounts for the 
uncertainty in the string tension a. 

In the third and fourth columns of Table 1, we show results from Refs.[21, 22] 
where anisotropic lattices have been used with different temporal and spatial spac- 
ings. The first error comes from the combined uncertainties from the continuum-limit 
extrapolation and from the anisotropy, the second from the uncertainty in a hadronic 
length scale playing a role similar to the string tension. 

It is worth pointing out that most of the investigations in lattice consider SU(N C ) 
theories at finite N c . However, in the AdS/CFT correspondence, the Yang-Mills 
theory is in the large N c limit. In order to compare the large- N c gauge theory to its 
finite counterpart, Ref. [23] calculated the lightest and the first excited scalar glueball 
masses in gauge theories for increasing N c and it was found a mass difference for the 
glueballs of about only 5% between the N c = 3 and the large- N c gauge theories (see 
also [26]). These results are shown in the last two columns of Table 1 where we used 
the mean value y/a = 440 MeV of [20]. 
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2.2 The Soft- Wall model in brief 

The AdS/QCD Soft- Wall model is defined by the action [12]: 

S = I d 5 if^ e-*£ , (2.1) 



where C is the Lagrangian density, $ is the static background dilaton field and g is 
the determinant of the metric tensor of AdS 5 , given by 

ds 2 = g mn dx m dx n = —{r 1 ^dx lx dx v + dz 2 ) , (2.2) 

with r/^ = (— , +, +, +) the Minkowski Ad metric. R is the AdS 5 radius and z is the 
holographic coordinate allowed to run from zero to infinity. In the following, we will 
set it! = 1 as it does not enter expressions of any four dimensional physical quantity 
as particle masses. In this model, the infrared cut-off is represented by a background 
dilaton field chosen as $ = cz 2 where the dilaton parameter c has the dimension of 
a squared mass. 

For the case of a massive scalar bulk field X = X(x, z), the action takes the form 

[13]: 

S = J d 5 x^-ge-^ [g mn d m Xd n X + m 2 AdS X 2 ] (2.3) 
and the equation of motion is 



dz{ ± e -H*) dzX ) + le"*( - ±e-*^m 2 AdS X = . (2.4) 

z° z° z° 

Representing the scalar field through a Ad Fourier transform: 

X(x,z) = J ^e**X(q,z), 



one has: 

d z (e~ B ^d z X) - q 2 e~ B ^X - ^fe'^I = (2.5) 



ml 



z 2 



with B(z) = <&(z) + 3\n(z). Under the change of function X = e B//2 V, one finally 
gets a Id Schrodinger-like equation: 

- d 2 z Y + V(z)Y = -q 2 Y (2.6) 

with the 5d effective potential: 

V(z) = ^-^ + ^ = c 2 z 2 + ^ + 2c+^f. (2.7) 

4 2 z 2 Az 2 z 2 

The normalizable solutions of Eq.(2.6) correspond to a discrete spectrum of Ad masses 
q 2 = — m 2 . The 5d mass is related to the scaling dimension A of the dual Ad boundary 
operator as 

m\ dS = A(A - 4) . (2.8) 
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The scalar glueball is associated with the local gauge-invariant QCD operator 
tf{G liV G tll/ ) defined on the boundary spacetime and which has classical dimension 
Adass. = 4 implying then rriAdS — 0. The corresponding mass spectrum was found 
in [13]: 

m 2 Gn =4c(n + 2) . (2.9) 

If one fixes c by the mass spectrum of the vector p mesons in AdS/QCD as c = 
0.2325 GeV 2 = (0.482 GeV) 2 [12], one obtains the results shown in Table 2. 



0++ 


0++* 




Q-| — |-*** 


1.364 


1.670 


1.929 


2.566 



Table 2: Scalar glueball masses in GeV from the Soft- Wall model [13]. 

One points out that these results from the Soft- Wall model are systematically 
smaller than the lattice values of glueball masses reviewed in the previous subsection. 
Furthermore, it is interesting to note that even changing the value of c, one can not 
fit the lattice glueball masses using the relation (2.9). 

We will see how these masses can be increased by the introduction of the anoma- 
lous contribution to the dimension of the scalar glueball operator. The anomalous 
dimension coming from quantum effects depends on the renormalization scale p of 
the four dimensional boundary theory. On the other side, holography implies that 
the fifth z coordinate representing the energy scale at which physical processes occur 
is inversely proportional to p [27]. So, the bulk mass (2.8) with A the full dimension 
will depend on the z coordinate, thus acting as an additional effective potential for 
the Schrodinger-like equation (2.6). 

2.3 The anomalous dimension of the scalar glueball operator 

We calculate the full dimension of the scalar glueball operator in our phenomeno- 
logical holographic model inspired by QCD. For this, following [28], we consider the 
trace anomaly of the QCD energy-momentum tensor: 

Tj^ = -^—-trG 2 + fermionic terms (2.10) 

^ 167TCT 

where the beta function is defined as usual as 

^ - si < 2 -"> 

with a = g 2 M /An and qym the Yang-Mills coupling constant. 

On the other hand, for any operator O, we have the following scaling behaviour: 

< 2 - 12 > 
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where the full dimension A c = A c ; ass . + 'y(fj) is given in terms of the classical 
dimension A c ; ass . and the anomalous dimension 'y(fj). In the following, the fermionic 
contribution on the r.h.s of (2.10) will not be considered, since we are interested 
only in the equation for the operator trG 2 . Thus, by taking into account the scalar 
glueball operator contribution of the QCD trace anomaly, Eq.(2.12) gives: 



8vra 2 



tr G' 



d 



din /j, \8na x 



tr G' 



(3(a) 



-((3'(a)--(3(a)-A G2 )^trG 



a ' ' ' %-na 2 



(2.13) 



where the prime denotes the derivative with respect to a. The trace scales 
classically, that means A t m = 4. This finally implies that the scalar glueball operator 
tr G 2 has the full dimension: 



A G2 = 4 + p'(a) - -p(a) . 

a 



(2.14) 



From now on, we will express the beta function in terms of the 't Hooft coupling 
A = N c g\ M = 4nN c a and the prime will now denote derivative with respect to A. 
The anomalous dimension defined above implies a correction for the bulk mass so 
that the potential (2.7) becomes 



15 



V(z) = c 2 z 2 + — + 2c+- 



4 + /3'(A)--/3(A) 



£'(A) - /(A) 



(2.15) 



In the following, we will solve the Id Schrodinger-like equation (2.6) with this 
generalized potential for different beta functions and determine the corresponding 
glueball mass spectra. 



3. Glueball masses from beta functions 

In this section, we will consider some beta functions discussed in the literature [29, 30] 
and adjust their different parameters in order - if possible - to obtain a mass spectrum 
compatible with the glueball lattice results. We will also demand that the various 
beta functions have a ultraviolet perturbative behaviour similar to QCD for small A 
at 1-loop order: 

/3(A) ^-M 2 , (3.1) 

where b is a positive constant in order to guarantee the asymptotic freedom property 
of QCD. 

Note that the r.h.s. of Eq.(3.1) was considered in [29] as an attempt for a beta 
function when investigating the heavy quark- ant iquark interaction potential in some 
Renormalization Group revised AdS/QCD models. However, in our phenomenologi- 
cal framework, such a beta function can not give additional anomalous contribution 
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to the effective potential, since the last term in (2.15) vanishes identically in this 
case. 

In the holographic models of QCD, the energy scale is related to the fifth coor- 
dinate of the AdS space by /i — 1/z. Hence, we can write an equation for the beta 
function in terms of z as 

^=/»<A(*0> =► *^ = -«AM). (3.2) 

We will choose the value of X(zq) = Ao at the energy scale given by the mass 
M z o of the neutral Z° boson so that 

z = l/M z o = 0.010966 GeV" 1 (3.3) 

and corresponding to the world average value of the strong coupling [31]: 

a s = 0.1184 A = 4.464 (N c = 3) . (3.4) 

3.1 Beta function with a cubic IR asymptotic behavior 

First, let us consider the following beta function [29]: 

f 3(X) = -b X 2 -b 1 X 3 (6o,6i>0). (3.5) 

This beta function behaves as —A 3 in the infrared. 

In this case, the renormalization group equation can be exactly solved, finding 

where X(z ) = A fixes the integration constant. 

On the one hand, when A 1, we recover, as it must be, a QCD-like behaviour: 

at the leading order. On the other hand, the limit A — > oo is reached for a finite 
value of z given by 

_ z e^ 

Zmax ~ bi/b 2 ■ V-O) 

This beta function implies that the potential of Eq.(2.15) takes the form: 

V(z) = c 2 z 2 + + 2c - ^6iA(^) 2 (4 - b^z) 2 ) . (3.9) 



Az 2 ■ 



So, in contrast with the Soft-Wall model, V(z) goes to infinity at z — z max . 
In other words, there is an infinite wall located at z max and the solutions of the 
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Schrodinger like equation will be non- vanishing only in the region < z < z max . 
As an unexpected consequence, the beta function (3.5) leads to the formation of a 
hard wall in the 5d effective potential. The fact that the background was initially of 
the Soft- Wall model does not have significant influence in getting the glueball mass 
spectrum in this shown in Table 3 where a vanishing dilaton parameter c is 

in fact a quite acceptable value. 

In order to do the numerical analysis for this beta function, we employ a standard 
shooting method and assume that the ratio between bo and b\ follows the pertubative 
relation bi = (51 /121)&q valid for pure SU(3) C gauge theory [18, 32]. So, we will con- 
sider c and bo as the independent parameters of the model. We use these parameters 
to fit the first two isotropic lattice glueball masses within their incertainty ranges 
[20] given in Table 1. Then, we calculate the other masses and compare them with 
the third and fourth masses listed in Table 1. It appears a lower bound for bo below 
which there is no value of c able to match the mass spectrum. When bo attains its 
lower bound b { min) = 0.03244, only the value of c = (0.385 GeV) 2 allows us to fit the 
masses shown in the first row of Table 3. 

Increasing the value of bo, we find a range of possible values for c. For example, 
for bo = 0.03250, we can fit the first two masses by considering values in the interval 
0.300 < y/c < 0.366 GeV. For b = 0.03270, we find the interval < y/c < 0.27 GeV. 
We also find that for c = 0, we can fit the first two masses provided that 0.03258 < 
bo < 0.03290. For the sake of illustration, in the last two rows of Table 3 we show the 
mass spectra for the mean values of y/c and bo corresponding to the ranges associated 
respectively with b = 0.03250 and c = 0. Finally, note that the values found for b 
are relatively close to the perturbative &[ l peri " )l . 



bo 




%max 


0++ 


++* 




Q++** 


0.03244 


0.385 


3.286 


1.570 


2.567 


3.545 


4.518 


0.03250 


0.330 


3.247 


1.546 


2.574 


3.570 


4.559 


0.03270 


0.000 


3.120 


1.515 


2.625 


3.677 


4.714 



Table 3: Masses for the scalar glueball groundstate and the excitations J PC = ++ for the 
beta function with a cubic IR asymptotic behavior for different values of bo and c. Masses 
and yfc are expressed in GeV. For completeness, we also show the values of z max expressed 
in GeV" 1 . 



3.2 Beta function with an IR fixed point at finite coupling 

The beta function: 

/3(A) = -bo\ 2 (l - £\ (V > 0) (3.10) 



^or instance, b { pert) ~ 0.03237 for n f = 5. 
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reproduces the perturbative /3(A) ~ —bo\ 2 at 1-loop order in the ultraviolet and 
behaves as /3(A) ~ +A 3 in the infrared. This beta function has been used for instance 
in Ref.[30]. Solving Eq.(3.2) for the beta function (3.10), we find: 

z(\) = Ce~^e^ ln{ \^ } (3.11) 

where C is an integration constant. Note that there are two solutions, one for A > A* 
and another for A < A*. In the first case, the solution does not show asymptotic 
freedom. So, we will consider only the physically relevant case A < A* for which we 
fix C imposing A(^ ) = Ao and we finally obtain: 

z( A) = 2 „e^-i)0^£)* (3.12) 

where the values of z Q and A are given in Eqs.(3.3) and (3.4). In particular, for 
A <C 1, we recover the QCD asymptotic behavior of Eq.(3.7). 

Writing A as a function of z and substituting into the beta function and then in 
the potential (2.15), we solve numerically the corresponding Schrodinger equation to 
find the glueball masses. The masses depend on the parameters: b , c and A*. For 
this case, we did not find a reasonable set of parameters that fits the masses for the 
isotropic lattice [20], the reason being essentially the gap between the groundstate 
and the first excited state (around 1.3 GeV, to be compared with the smaller gap of 
0.94 GeV for anisotropic lattice [21]). However, for the anisoptropic lattice results 
[21], we can fit the first two masses with the choice &o — 0.091, y/c = 0.55 GeV and 
A* = 600. These results and also the predictions for the third and fourth excited 
states are shown in Table 4. It is interesting to remark that we can still fit the first 
two masses keeping two of these parameters with the values shown above and varying 
the third one in the intervals: 0.03902 < b < 0.04005, 0.535 < y/c < 0.616, GeV 
and A* > 552. The ranges of values for b and y/c turn out to be rather stringent. As 
for the lower bound on A*, it allows the IR fixed point to be pushed away to infinity, 
giving rise then to a beta function with a fixed point at infinite coupling. Such beta 
functions have been considered, for instance, in [17]. 



bo 




A* 


0++ 






Q-| — |-*** 


0.0391 


0.55 


600 


1.734 


2.380 


3.053 


3.727 



Table 4: Masses for the scalar glueball states for the beta function with an IR fixed point 
at finite coupling. Masses and y/c expressed in GeV. 

3.3 Beta function with a linear IR asymptotic behavior 

Now we consider the following beta function [29, 33]: 

f}(\) = -J£L. (&oA>0) (3.13) 
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that behaves like the perturbative QCD beta function at 1-loop order and decreases 
asymptotically as —A in the infrared. 

In this case, the 5d effective potential reads: 

V ^~ CZ + 4z 2 + ZC+ Z 2 (l + M^)) 4 ^ J 

and the solution of the renormalization group equation is 

X(z) = y—^ (3.15) 



biW ^te (f } Wh j 



where is the Lambert function (satisfying the property W(x)e w ^ = x) and 

the integration constant is fixed such that X(z ) = A . The potential can then be 
put into the form: 



15 lb K i+w( ^ ( ^ )6o/bi) ) + ^ 

v{z) = cV + if + 2c + ^ 1A T S= V~ ' ; :u { ] 

1 [i + ^(^(?) 6o/bl ) 

In the infrared limit (when z is large), the coupling constant behaves as A(z)| 2>20 
z b °/ bl such that 



^)U» 20 ~c z + — + 2c+ ^ ^ ■ 



In other words, when z — )■ oo, the potential contribution stemming from the anoma- 
lous dimension is subleading with respect to the oscillator-like term c 2 z 2 coming form 
the soft wall background dilaton field. Consequently, the Regge-like behaviour of the 
mass spectrum m 2 ~ n for large n is not modified, in contrast to the case of the beta 
function with a cubic IR asymptotic behavior (cf. subsection 3.1) which led instead 
to a hard wall-like mass spectrum m 2 | n>i ~ n 2 . 

The masses depend on the parameters: b , b\ and c. Once again, we have 
not been able to find a reasonable set of parameters that fits the masses for the 
isotropic lattice but, as for the anisotropic lattice results, we can fit the first two 
masses with the choice bo = 0.0395, y/c = 0.49 GeV and b\ = 0.001. These results 
and also the predictions for the third and fourth states are shown in Table 5. The 
first two masses can still be fitted keeping two of these parameters with the values 
shown above while varying the third one within the intervals: 0.0391 < bo < 0.0404, 
0.396 < y/c < 0.593 GeV, and < h < 0.00156. When 6 X goes to zero, the potential 
presents the form of a very deep and sharp well mimicking a hard wall-like potential, 
at least for the lowest-excited states. 
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61 


0++ 


Q++* 


Q++** 


Q-| — |-*** 


0.0395 


0.49 


0.001 


1.704 


2.485 


3.285 


4.096 



Table 5: Glueball masses in GeV for the the beta function (3.13). \fc is expressed in GeV 
whereas bo and 61 are dimensionless. 



4. Effective beta functions from glueball masses 

The glueball masses obtained from lattice studies shown in Table 1 lie approximately 
in linear Regge trajectories, namely, the square of the masses has approximately a 
linear relation with the radial quantum number. In this section, we build up beta 
functions that reproduce this behaviour. In particular, we will consider the masses 
from the isotropic lattice calculations shown in the second column of Table 1. A 
linear fit for these masses corrresponds to 

m 2 n = 4.50n + 2.51, (4.1) 

where the masses squared are given in GeV 2 . This linear fit implies the glueball 
masses given in Table 6. The first mass is just 1% above the lattice result induing 
the error bar. The other three masses are within the error bars. 



0++ 






Q-| — |-*** 


1.585 


2.648 


3.393 


4.001 



Table 6: Glueball masses in GeV from the linear fit of Eq.(4.1). 



In order to obtain the class of beta functions that lead to this Regge linear 
spectrum, we consider the soft-wall potential modified by the contributions from the 
anomalous dimension of the scalar glueball operator (2.15): 

V(z) = c 2 z 2 + il + 2c + 1 [4 + /(*)] /(*) (4.2) 



4z 2 z 2 

where 



It is convenient to rewrite f(z) in terms of the coupling X(z): 



x A' A 

-1-4 (in (4-4) 



dz V A 2 
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where prime denotes derivative with respect to z. 
We use the following ansatz 



{f{z) + 2f = k + k lZ 2 + k 2 z\ (4.5) 
where k, k\ and k 2 are constants. The above ansatz requires the positivity condition 

k\ < Akk 2 . (4.6) 

The potential takes the form 

k - - 

V(z) = (c 2 + k 2 ) z 2 + + (2c + h) (4.7) 

Solving the Schrodinger equation (2.6) with this potential for k = m 2 where m is an 
integer, one finds the 4d mass spectrum: 



A^/c 2 + k 2 n + + k 2 (m + 1) + 2c + k x . (4.8) 



This equation can be matched to the glueball linear fit of Eq.(4.1). This matching 
and the positivity condition (4.6) determine the allowed ranges of the parameters ki, 
k 2 , c, and m. It turns out that c m i n < c < c max where 

Cmin = -(0.983 GeV) 2 ; c max = (0.344 GeV) 2 . (4.9) 

The extreme values of c saturate the positivity condition, i.e., k\ = 4kk 2 . We also 
observe that k\ takes only negative values. 

Let us now determine the corresponding beta functions. From Eqs.(4.4) and 
(4.5), we have to solve the differential equation: 



d_ 

dz 



In (^|) =- z T- z Vk + k lZ 2 + k 2 z* (4.10) 



To reproduce QCD asymptotic behaviour at leading order shown in Eq.(3.7), we 
impose that 

lim (-L-Y = , (4.11) 

which implies: 

i. d , fX'(z)\ 1 

Jas 1 "^]"-;- (412) 

Then we must have k = 4 (and m = 2). For small z, we also must take the upper 
sign in Eq.(4.10). Note that this analysis of the perturbative small z limit does not 
impose the sign of Eq.(4.10) for large values of z as it will be discussed below. 

We will treat separately the cases k\ = 16k 2 and k\ < lQk 2 , since then the 
differential equation (4.10) assumes different forms. 
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4.1 Case k\ = I6k 2 

In this case, c can take any one of the two values c m i n and c max - Then, the ansatz 
(4.5) and the differential equation (4.10) together with (4.4) reduces to: 

Since k\ is negative, there is a critical value of z given by z\ = —8/ki, such that the 
r.h.s. of Eq.(4.13) is zero. We get z c = 1.8615 GeV" 1 and z c = 1.337 GeV" 1 for 
c = c m in and c = c max , respectively. For z < z c , the sign of (4.13) is fixed by the 
perturbative limit. For z > z c we are free to choose the sign, without spoiling the 
smoothness of the coupling X(z). In the following, we are going to consider the two 
possibilities. 

4.1.1 Solution without change of sign at z c 

We take the upper sign of (4.13) for all values of z and find an analytical solution: 

x{z) = i-b ^[m-^)-E i (-^4)} (414) 

where Ei{x) is the exponential integral function and the initial condition X{zq) = A 
has been used. We assume from now on that b takes its perturbative expression 
with iV c = 3: 

8tt 2 V 3 9 J K J 

This equation introduces the number of flavours rif. However, the four lattice glueball 
masses that we use in this section were obtained from SU(3) C with no flavour. Note 
that unquenched lattice QCD only provides the ground-state scalar glueball mass, 
usually suffering from severe computational difficulties (lack of high statistics, coarse 
lattice spacing, etc) [34, 35, 36, 37]. By varying rif in the coupling (4.14) we will 
describe effectivelly the infrared contribution in b appearing in (4.11). We will see 
in the following that the low energy behaviour of the coupling changes according to 
n f . 

The beta function for this particular case reads: 

ftz) = -zX'(z) = — . (4.16) 

\l-b^[Ei{-^)-Ei{- k -±zl)]} 

Considering the two possible values for c which fixes k\ and k 2 , we have plotted in 
Figure 1 the couplings and the beta functions in terms of the holographic coordinate 
z for 71/ = 0, 3 and 6. Both X(z) and j3(z) diverge at a finite value of z — z max where 
the denominators in (4.14) and (4.16) vanish. Note that z, max depends on n/. In 
Figure 2, we show the dependence of the beta functions on the coupling A for these 
two values of c. 
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0.0 0.5 1.0 1.5 2.0 2.5 3.0 -200 



Figure 1: The 't Hooft couplings X(z) and the effective beta functions f3{z) as a function 
of z as defined in Eqs.(4.14) and (4.16) for different values of nj and c = c max (upper 
panels) and c = c m i n (lower panels). 




Figure 2: The effective beta function /3(A) in terms of A for different values of nj when 
c = c max (left) and c = c min (right). 



4.1.2 Solution with change of sign at z c 



Now we choose for the differential equation (4.13) the following form 



d ( X' 



-(1 + ^ 2 ) 
(3 + *£) 



z < z c , 
z>z c . 



(4.17) 
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The solution is 



\(z) = (l-e{z-z c ))X < {z) + 



e(z-z c ) 



(4.18) 



A<(z c ) 



1 



+ |(| + l)e : 




b e 



yielding the beta function: 



p( z ) = (i-e(z-z e ))p < (z) 



1 



9(2 - z c ) 



2- 




6 e 




(4.19) 



A<(^ c ) 



where O(z) is the step function and \<(z) and < (z) are the analytical solutions 
(4.14) and (4.16) of the previous subsection. 

We observe that the shapes of the 't Hooft coupling and the beta function are 
similar to the preceeding case, shown in Figure 1, diverging at z — z max . 

Depending on nj, z max can be smaller or larger than z c . When z max < z c , only 
the solution \<(z) is present. This happens when n/ < 2 (negative c) and n/ < 1 
(positive c). When z max > z c , we have another solution, given by the whole r.h.s. of 
(4.18), and which gives rise to a larger value of z max with respect to the preceding 
case. 

Below, we display in Figure 3 the dependence on A of the beta function for three 
values ofnf and the two possible values of c. 



Figure 3: The beta functions /3(A) as a function of A for different values of rif and the 
two possible values c ma x (left) and c m i n (right). 

Let us emphasize that changing the sign in (4.13) leads to two different solutions. 
Nevertheless, the data we have used do not allow us to select a unique beta function. 

4.2 Case kj < 16k 2 

In this case, the square root on the r.h.s. of (4.10) does not vanish for any value of 
z. Therefore, the upper sign is the only possibility here. Then, we find out: 



1 V = b ( V4 + fcg + gg + ^ 2 + 2) (fcx + AVh) i J , +klZ , +k2Z , 






(4.20) 
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which cannot be solved analytically. The parameter c calculated from Eqs.(4.1) and 
(4.8) with m = 2 can assume any value inside the range c m i n < c < c max . 

In order to analyse the behaviour of the general equation (4.20), let's write its 
integral from z to z in the form: 



1 

Then, the 't Hooft coupling reads 

\{z) 



1 



X(z) 



F(z,z ) 



1 - X F(z,z ) 



(4.21) 



(4.22) 



For small enough values of n/, namely rif < 4, there is always a z = z max at 
which A (z) diverges. Moreover, this occurs for any value of c. Then, the shapes of 
X(z), f3(z) and (3(X) are similar to the ones displayed in the Figures 1, 2 and 3. 

On the contrary, for n/ > 5, we will see that the results turn out to be completely 
different. In the sequel, we will consider the illustrative case rif — 6. 

First, we observe that for c close to the extreme values of its allowed range, the 
coupling solution remains singular. Like the two cases considered above, there is a 
value z = z m ax for which the denominator 1 — XoF(z max , z ) vanishes. 

However, there is an interval for c, namely 



c e [-(0.957 GeV) 2 , (0.164GeV) 2 



(4.23) 



for which this denominator never vanishes. As a result, the coupling and the beta 
function are finite for all z. We show in Figure 4 the behaviour of 1 — XqF(z, zq) = 
Xq/X(z) as a function of z and c. 




Figure 4: The function Xq/X(z) and the zero-value plane surface. The coordinate z and 
the parameter c run respectively from to 50 cind. from its extremal values c m i n — -0.966 
GeV 2 to c max ~ 0.118 GeV 2 . 



Further, in Figure 5, we show the 't Hooft couplings and the beta functions for 
some noticeable values of c. The red and green curves correspond to the two limiting 
values c = (0.165 GeV) 2 and c = —(0.958 GeV) 2 beyond which X(z) and j3(z) are 



singular at z max . On the other hand, for any value of the parameter c inside this 
interval, the 't Hooft coupling presents an IR fixed point A*. Accordingly, the beta 
function vanishes for sufficiently large z. Figure 5 also displays the case c = (no 
background dilaton field) and the mid-value c = —(0.687 GeV) 2 as well. 



A(z,c) filf) 




Figure 5: The 't Hooft coupling X(z) and the associated beta function f3(z) in terms of z 
for some remarkable values of c in the case nj = 6. 



c (GeV 2 ) 


Zmax (GeV) 


A* 


(0.165) 2 


4.300 




(0.164) 2 




61670 







300 


-(0.686) 2 




38 


-(0.957) 2 




3005 


-(0.958) 2 


4.873 





Table 7: Numerical estimates for z max and A* for some values of the dilaton parameter 
inside the interval (4.23) and immediately beyond. 

In Table 7, we list the numerical estimates for z max and A* for some characteristic 
values of c. Note that A* takes very different values. 

Finally, we plot the beta functions in terms of A in the Figure 6. Within the 
interval (4.23), /3(A) shows the typical behaviour associated with an IR fixed point 
at finite coupling. On the other hand, for c close to its extremal values, we recover 
the monotonous decreasing of /9(A). 
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Figure 6: beta functions in terms of A for rif = 6 and different values of c when kf < 16^. 

5. Conclusion 

The AdS/CFT correspondence provides a well-defined relation between the dimen- 
sions of Yang-Mills theory operators and the masses of the corresponding supergrav- 
ity fields. Nevertheless, for QCD-like gauge theories, the dimensions of operators 
receive anomalous contributions coming from quantum effects. In this paper, we 
have developped an extended AdS/QCD approach taking into account the anoma- 
lous dimensions of the operators entering the expression of QCD trace anomaly. We 
considered the scalar glueball spectroscopy within the AdS/QCD Soft- Wall frame- 
work. Following this extended model, the mass of the scalar bulk field acquires a 
dependence on the gauge theory beta function. This modification is then reflected 
in the eigenvalue equation that determines the 4d glueball masses. 

We investigated the relation between beta functions and masses in two differ- 
ent (complementary) ways. In Chapter 3, we considered some possible QCD beta 
functions discussed in the literature and calculated the corresponding glueball mass 
spectra, comparing with lattice results. For each case, we found ranges of parameters 
that fit lattice data. Interestingly, some beta functions give rise to the formation of 
a hard-wall in the 5d effective potential. 

In Chapter 4, we followed the inverse way: we started from the glueball masses 
calculated using lattice and found out beta functions that lead to this spectrum. We 
find different types of beta functions classified according to their IR behaviour. For 
a class of beta functions, an IR cutoff, z max , which is not present in the Soft- Wall 
model naturally emerges and gives rise to asymptotically decreasing /3(A). There is 
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also another class of beta functions that vanish at large A, leading to an IR fixed point 
at finite coupling. The IR vanishing beta functions offer a large range of possible 
values for the IR fixed coupling. 

A better knowledge of mass spectroscopy in different channels would further 
constrain the holographic parameters. This could help answering the question of the 
presence or not of an IR fixed point and its value. 
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